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ABSTRACT 

We consider the simplest instabilities involving multiple unstable electrostatic 
plasma waves corresponding to four-dimensional systems of mode amplitude 
equations. In each case the coupled amplitude equations are derived up to third 
order terms. The nonlinear coefficients are singular in the limit in which the 
linear growth rates vanish together. These singularities are analyzed using tech- 
niques developed in previous studies of a single unstable wave. In addition to the 
singularities familiar from the one mode problem, there are new singularities in 
coefficients coupling the modes. The new singularities are most severe when the 
two waves have the same linear phase velocity and satisfy the spatial resonance 
condition /c2 = 2fci. As a result the short wave mode saturates at a dramatically 
smaller amplitude than that predicted for the weak growth rate regime on the 
basis of single mode theory. In contrast the long wave mode retains the single 
mode scaling. If these resonance conditions are not satisfied both modes retain 
their single mode scaling and saturate at comparable amplitudes. 
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1 Introduction 



Recently we presented a detailed analysis of the amplitude equation for a single unstable 
electrostatic mode in an unmagnetized Vlasov plasma (henceforth (I)).||l| The analysis reveals 
a fundamental difficulty with the derivation of amplitude equations for this class of problems: 
the coefficients in the amplitude equations become singular in the limit in which the growth 
rate 7 of the unstable wave is allowed to vanish. Although these singularities can be removed 
by an appropriate 7-dependent rescaling (see below) the analysis shows that amplitude 
equations of this type cannot be truncated at any finite order. [0, R[ Nonetheless the scaling 



identified by the theory predicts the amphtude at which weakly growing waves will saturate, 
and hence is of fundamental importance both in plasma physics and in the closely related 
problem of shear flow instability of ideal fluids. 

In view of the importance of the predicted scaling for applications we investigate here the 
effects of including additional unstable modes. We consider only the simplest possibilities, 
those requiring a four-dimensional system of amplitude equations. There are three such 
instabilities distinguished by the symmetry of the equilibrium and whether the unstable 
modes have real or complex eigenvalues. Wc find that unless the two modes satisfy a strong 
resonance condition the presence of the second mode does not alter the saturation amplitude 
of the original mode. However, in the important resonant case in which the phase velocities 
of the two modes are the same and their wavenumbers ki, /c2 satisfy k2 = 2ki the long wave 
mode suppresses dramatically the saturation amplitude of the short wave mode. 

The amplitude equation defines, in the limit 7 — > 0"*", a kind of singular perturbation 
problem whose detailed features reveal asymptotic scaling behavior of the nonlinear wave. 
This is a key idea behind our approach and the reader is referred to (I) for a more detailed 
discussion. For the single mode instabilities there is only one amphtude A{t) and one seeks 
a scaling A(t) = 'j^a^'jt) such that the evolution equation for a(r), r = 7^, has a non- 
singular limit as 7 ^ 0+. In dissipative problems, the critical eigenvalues are isolated on 
the imaginary axis, and /3 = 1/2 is the generically expected exponent. As a result, in the 
generic case, the amplitude equation can be truncated at third order. This is not so for an 
unstable electrostatic wave. Here the situation is quite different because the Vlasov equation 
is Hamiltonian and the eigenvalues of the mode merge with a continuous spectrum on the 
imaginary axis at criticality, i.e. as 7 ^ 0"*". As a consequence the nonlinear coefficients 
in the amplitude equation are singular as 7 ^ 0"*" signalling strong nonlinear effects that 
saturate the unstable linear growth at exceptionally small wave amplitudes. A quantitative 
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signature of this reduction of the nonhnear wave amphtude is a larger exponent: /? = 5/2 
for plasmas with multiple mobile species and /3 = 2 in the limiting case of infinitely massive 
(fixed) ions and mobile electrons. In fact, the analysis presented in (I) showed that setting 
/9 = 5/2 yielded a theory that was finite to all orders in the amplitude expansion as 7 ^ 0^. 

In this paper we investigate the coupled amplitude equations for two unstable modes 
with amplitudes A{t) and B{t). The coupled equations contain the single mode instabilities 
due to excitation of only A or only B, and the previously studied singular coefficients govern 
these special cases. There are now separate scalings possible for each amplitude, 

|A(t)|=7^a(7it) |S(t)|=7^^&(72t), (1) 

and we know (3j > 5/2, j = 1,2 is required to control the singularities in the single mode 
coefficients. We seek to determine the possible singularities in the coupling coefficients 
between A and B that are new; specifically we wish to know if these singularities can 
dominate the single mode singularities and require new nonlinear scalings for the instability 
with two simultaneously growing waves. 

In the remainder of this introduction, we summarize our notation and state some relevant 
results about the linearized theory. Section || enumerates the possible instabilities described 
by four-dimensional systems of amplitude equations. In each case the general form of the 



amplitude equations can be anticipated on the basis of symmetry. Section p.3| describes 
the procedure for calculating the coefficients in these amplitude equations and summarizes 
the results for the leading terms through third order. The singularities that arise in these 
terms in the limit of weak instability are analyzed in Section ^, and their consequences are 
discussed in Section 4. The paper ends with a brief conclusion. 
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1.1 Notation 

We briefly summarize the notation of (I) which will be used. Our model is the one- 
dimensional, multi-species Vlasov plasma defined by 



Qpis) gf?{s) gp{s) gj? .00 



Here x, t and v are measured in units of w/cje, uj~^ and m, respectively, where m is a cho- 
sen velocity scale and uol = ATre^ne/nie. The plasma length is L with periodic boundary 
conditions and 

dx r dvF^'\x,v,t) = ( ^) L, (3) 

J-L/2 J-00 \ Tie J 

where Qs = ezg is the charge of species s and = qsnie/enis. Note that k*-^-* = —1 for 
electrons and that the normalization (j^) for negative species makes the distribution function 
negative. 

Let Fq{v) and f{x,v,t) denote the multi-component fields for the equilibrium and per- 
turbation, respectively, and k the matrix of mass ratios. 



Fo = 



V 



V J 

then the system (0) can be concisely expressed as 



K = 



( ■■■ \ 

■■■ 

V : : : ; 



(4) 



(5) 



where 



l=—oo 

/ = 

vfiiv) + K ■ r^iiv) Es' I-oo dv' ft'\v')] I ^ 0, 



-il 



(6) 
(7) 
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with rii{v) = —d^Fo/l"^, and 

^ f .. ^frr 

m=—oo ;=— oo 

In the spatial Fourier expansion @, / denotes an integer muhiple of the basic wavenumber 
27t/L, and a primed summation as in (H) omits the / = term. The notation k ■ rii{v) or 
■ d^fm-i denotes matrix multiphcation. 
An inner product is needed in Section ^ to derive the amphtude equation. For two multi- 
component fields of (x, v), e.g. B = {B^'-'\ B^'^\ B^'^\ ..) and D = {D^'^\ D^'^\ D^'^\ ..), we 
define their inner product by 



l-L/2 rOD /-L/2 

{B,D) = dx dvB^'\x,vYD^'\x,v)= dx<B,D>, (9) 

, J-L/2 J-oo J-L/2 



where 



/CO 
dvB^'\x,vyD^'\x,v). (10) 

1.2 Summary of linear theory 

The spectral theory for C is well established, and the needed results are simply recalled to 
establish our notation. ^, |^ The eigenvalues A = —ilz of C are determined by the roots 
Ai{z) = of the "spectral function'' 



A,w^i+rrf.^-^'"'"'"<^">. (11) 

J-oo V — Z 

The linear dielectric function ei{z) is obtained on replacing the contour in (0) by the Landau 
contour for lm{z) < 0; for lm{z) > 0, Ai{z) and £i{z) are the same function. 

Associated with an eigenvalue A = —ilz is the multi-component eigenfunction \E'(x, t;) = 

^iix ip[y'^^ where 

i^{v) = (12) 



There is also an associated adjoint eigenfunction ^(x, f) = e^^^%lj{v) / L satisfying (\^, \E') = 1 
with 

^(^) = -TTTzJ——-^- (13) 



AJ(z)*(t; - z*)' 

Note that all components of il){v) are the same. The normalization in ([T3|) assumes that 
the root of A/(z) is simple and is chosen so that Kip^ip >= 1. The adjoint determines the 
projection of f{x,v,t) onto the eigenvector, and this projection defines the time-dependent 
amplitude of ^f, i.e. A{t) = (^,/). 

In section some of our results are conveniently stated in terms of the resolvent oper- 
ator, Ri{w) = {w — Li)~^ , whose general form follows from (J^) by solving [w — Li)f = g for 
/.[H, 1^ Here both g{v) = {g^^'^'' (v) , g^^'^\v) , . . .) and / are multi-component fields, and Ri{w) 
acts by 



Ri{w)g 



where 



il{v — iw/l) 



( (Riiw) gY^^Kv) \ 
(Riiw) gY'^^v) 



V 



Ai{lW/l) ^^-oo 



(14) 



dv 



v' — iw/l 



(15) 



2 Amplitude equations: general features 



Each of our instabilities can be formulated within a general framework as follows. The 
wavenumbers of the unstable modes '^2{XiV) and '^i{x,v) are given by \k2\ > \ki\ > 0, and 
the corresponding eigenvalues are Xj = —ikjZj for j = 1, 2, where Ak^{zj) = 0. With periodic 
boundary conditions, each wavenumber is a multiple of the minimum k, i.e. kj = 27mj/L 
with integer nj. We assume all roots are simple, i.e. A'f^,{zj) ^ 0; in addition, in the limit 
lm{z) — 5> of weak growth rates, the equation for the root is given by 



j:^^'H'iv,)=o, i+p/ 



dv 



0, 



V — Vi 



(16) 



where vj = R.e{zj) is the phase velocity at criticahty. These relations will be important for 
our analysis of the singularities in the nonlinear coefficients. 

The four-dimensional eigenspace is spanned by "^l, \l/2, ^2}' ^^"^ components 
of the distribution function along the unstable eigenvectors are identified by writing 



f{x, V, t) = v) + 5(t)^2(a;, v) + cc] + S{x, v, t), 



(17) 



where A{t) = /) and B{t) = (^^2, /) are the mode amplitudes and (^'j, S) = 0. In (|T 
= exp{ikjx) ipjl L is the appropriate adjoint function for Zj from ([T3|) . 
In the (A, B, S) variables, the Vlasov equation (|^) becomes: 



A 



B 



X2B + {^2,W)) 



(99 r ~ 

— = cs + ^^{f)-[{^^,^^{f))^, + {^2,^^{f))^2 + cc 



where 



l=—oo s' 



dv'fi'\v') 



(18) 

(19) 
(20) 

(21) 



In writing (0) we have used the adjoint relationship j, CS) = {C'^'^j, S) = Xji'^j, S") = 
and in ( pT] ) an integration by parts shifts the velocity derivative onto ip. 

These coupled equations are equivalent to (|^); however, by restricting them to the unsta- 
ble manifold we obtain autonomous equations for A{t) and B{t). The details of this restric- 
tion are discussed in an earlier paper and also in (I). The unstable manifold is tangent to 
E'^ at the equilibrium, and near Fq it can be described by a function H{x, v, A, A*, B, B*). 
Thus 



r(x, V, t) = [A(t)^i(x, v) + B{t)^2{x, v) + cc] + H{x, V, A{t), A*{t), B{t), B*{t)) (22) 
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represents a distribution function on W^. In this expression, the evolution of S is determined 
from H, i.e. 

/ H^''\x,v,A,A*,B,B*) \ 



S{x, V, t) = H{x, V, A{t), A*{t), B{t), B*{t)) 



H^'^Hx,v,A,A\B,B* 



\ 



(23) 



When this representation is substituted into (|T^) - (pOf ) we obtain 

A = \,A + {^^,^^in) 

B = \2B+i^2,W)) 



dS 
'dt 



r 



(24) 
(25) 
(26) 



Note that equations (|2^) - ( pSf ) define an autonomous flow describing the self-consistent 
nonlinear evolution of the unstable modes; this is the four-dimensional system we study. 

Symmetries of the Vlasov-Poisson system (^ and the equilibrium Fq{v) are important 
qualitative features of the problem. Spatial translation, 7^ : {x,v) {x + a,v), and re- 
flection, 71 : {x,v) — > (—X, —v), act as operators on f{x,v,t) in the usual way: if a denotes 
an arbitrary transformation then {a f){x,v) = f{a^^ {x,v)), where {a f){x,v) denotes the 
transformed distribution function. The operators C and A/" commute with 7^ due to the 
spatial homogeneity of Fq, and if Fq{v) = Fq{—v), then C and Af also commute with the 
reflection operator TZ. With periodic boundary conditions, x is a periodic coordinate so Ta 
and TZ generate 0(2), the symmetry group of the circle. If only the translation symmetry is 
present, then the group is S0(2). 

The action of translation Ta on f{x,v,t) implies an action on the variables {A,B,S). 



From ([T7|) we note that Taf{x, v, t) = f{x — a, v, t) is equivalent to 



Ta (A, B, Six, v)) = (e-^'^^M, e-'^^'^B, S(x - a, v)). 



ik2 



(27) 



The representation of TZ in the variables [A, B, S) depends on specific details of \I'i and 



in the individual cases discussed below. These symmetries determine the general form of the 
amplitude equations ( p^ - (^) for each of the instabilities we consider. 

2.1 Instability without reflection symmetry: Fq{v) ^ Fq{—v) 

When Fo{v) lacks reflection symmetry, the generic four-dimensional problem arises for modes 
with unequal wavenumbers k2 > ki > and complex eigenvalues. The roots, Ak-{zj) = 0, 
determine the phase velocities Vj = ujj/kj and the growth rates 7^ of the linear modes from 
the real and imaginary parts of the eigenvalue Xj = —ikjZj = 7^ — iuj. The corresponding 
eigenvectors are 

^j{x,v) = e'^^^'^lj^iv) = e''^' (-^ir^) > J = 1.2. (28) 

In this case, the identities A^^. (zj) = A_k.{zj) and Ak.{zj)* = A^^, (2;*) imply three addi- 
tional eigenvectors: \E'*, $j and $*, where 

^j{x,v) = e^'^^^jiv)*. (29) 

These eigenf unctions correspond to eigenvalues A*, —A* and — Aj, respectively, and fill out 
the eigenvalue quartets characteristic of Hamiltonian systems. In the absence of refiection 
symmetry, the eigenvalues are typically simple and the four-dimensional unstable subspace is 
spanned by {^^i, "^l, \E'2, ^2}- ^ beam-plasma system with a weak beam is the prototypical 
example of this instability. 

Since Fo{v) is spatially homogeneous, the amplitude equations ( ^4]) - ( pSj) always have 
translation symmetry (|2^) and we can apply standard results on the form of such symmetric 
equations. |0, ^ Hence we know the right hand side of (|^) - (|25D takes the general form, 

A = r{a)A + s{a){A*)''^-^B''' (30) 
B = p(a)5 + g(a) A"2(5*)"i^\ (31) 
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where a = (cxi, ct2, cts, (T4) denotes the four basic invariants 

ai=|Ap a2 = |5|2 (T3 = (T4 = (32) 

The integers rij refer to the wavenumbers kj = 27inj/L.^ The complex- valued functions 



r, s,p, and g are determined from (^If) - (|25D , but they can depend on the amplitudes only 



through the four invariants a"i, (T2, as and (J4. 

2.2 Instability with reflection symmetry: Fq{v) = Fo{—v) 

When the equilibrium is reflection-symmetric the transformation TZf{x,v,t) = /(— x, —v,t) 
commutes with C and A/" in the Vlasov equation. In this setting, the roots of A; (2;) can be 
either imaginary or complex, depending on the detailed form of Fo{v), and correspondingly 
we may encounter instabilities due to either real or complex eigenvalues. 

2.2.1 Real eigenvalues 

The description for real eigenvalues is quite similar to the examples without symmetry and we 
use the same notation (|28| ) for the unstable eigenvectors \E'i(a;, f) and 'if2{x,v). In this case, 
however, both linear phase velocities are zero. In addition, each of the real eigenvalues Ai 
and A2 has multiplicity two since the states TZ'^i = \E'^ and TZ'$2 = ^2 ^Iso eigenvectors. 
The translation symmetry (^) still holds, as well as the reflection transformation given by 

7^ {A, B, S{x, v)) = {A\ B\ S{-x, -v)). (33) 



The form of the amplitude equations (^) - ( pSj ) is the same except that the functions 
r,s,p and q in ( PDD - (^TJ) are now real-valued and depend on only three invariants Ui, (J2, 
and (T+ = (T3 + (74. This model applies for example to a reflection-symmetric two-stream 
distribution with instability at two wavenumbers. 



10 



2.2.2 Complex eigenvalues 

An instability in a reflection-symmetric system with complex eigenvalues likewise yields 
eigenvalues of multiplicity two. If only one wavelength is involved these result in a four- 
dimensional problem, and provide an additional example of what may be called a one-mode 
instability, cf. (1). Such an instability arises, for example, in a beam-plasma system with 
counter-propagating beams. [|1^] 

In the notation of our general framework, we let fci = — A;2 = ^ > and Zi = —Z2 = Zq; 
then Ak{zo) = implies A„fc(— zq) = 0, and these roots correspond to reflection-related 
eigenvectors for the eigenvalue A = —ikzo, i.e. A = —ikiZi = —ik2Z2: 

V v-zq ) 

^2{x.v) = {n^,)(x,v) = e-''''M-v) = e-''''(-^^^]. (35) 

V V + Zo J 

These solutions describe oppositely propagating linear waves with phase velocities vi = 
—V2 = uj/k, where zq = Vi + ij/k. In contrast to the previous examples, where both zi and 
Z2 sit in the upper-half plane (corresponding to positive wavenumbers) , in this case Z2 falls 
in the lower-half plane. 

Now the mode amplitudes in (|T7t) transform according to (^71) and applying TZ to ( |T7D 
yields 

n{A,B,S{x,v)) = {B,A,S{-x,-v)). (36) 
The 0(2) symmetry generated by TZ and Ta implies amplitude equations of the form 

( ^ ) = ^(/^i'^2,/i2) (yB^ +Q(/^i,/^2,/i2) ^* ) ' (37) 
where P and Q are functions of the invariants fii = \A\'^+\B\'^ and /i2 = AB with P(0, 0, 0) = 
A and Q{0, 0, 0) = 0. In this case a further simplification is possible; terms that do not 
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commute with the "phase-shift" symmetry, {A,B) {e^'^A,e^'^B), can be removed by near- 
identity coordinate changes, {A, B) = {A' + ^i{A', B'), i?' + $2(A', B')), to obtain the normal 
form: 



A 
B 



^(/^i,yU3) 



A 
B 



+ S{fi„fis)[\A\'-\B\'] 



21 / A 

-B 



(38) 



where /is = — l-Bp)^ and we have dropped the primes on H 
2.3 Amplitude expansions 



We wish to study the nonhnear terms (|2l|) in the amphtude equations - (^51) . The 
Fourier components /; follow from (^) 

fnv) = AMv)\k, + A*M^y6L-k, + B^2{v)6Lk, + B*M^yk-k2 + Hi{v), (39) 
and the amplitude expansion of Hi{v) begins with second-order terms, 

Hi{v) = [h{v)\A\^ + h2{v)\B\^]6i,o + h{v)BA*6i,k,-k, + h{vrB*A6i,k,.k, (40) 
+h4{v)A^6i,2k, + h4{vy{A*)^6i,^2k, + h{v)AB6i,k,+k, + h{vyA*B*6i,-k2~k, 
+h{v)B^6i,2k, + h{vy{B*y6i,-2k, + 0{3). 



Thus the nonlinear terms (21) can be written out in terms of the coefficients hi{v) in 
neglecting terms that are higher than third order in the amplitudes: 

{^l,^^{n) = riiO)A\A\^ + r2(0)A|5|2 + 4,,2fc,s(0)A*S + 5k,,3kA0)iA*YB (41) 
+6k„-k, [^2(0)^^5 + (gi(0) + P3(0))5*|A|2 ^ Q3(o)A*(5*)2 ^ g^(o)5*|fi|2l 

with Taylor coefficients 



ri(0) 



< dyijji, K- {hi - hi) > +— < d^tfji, K-ipl > 



(42) 
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r2(0) 



< dyi^i, K-h2> , < K-{hl-h5)> r* ~ 



ki + k2 



k2 — ki 



(43) 



.(0) 



P2(0) 

Qi(0) + P3(0) 



—I 



< dyilJi,K ■ ipl >/k2- < dyipi,K ■ iIj2 >/ki if k2 = 2ki 



—I < d^ipi, K, ■ hi > /k2 — < d^ipi, K ■ hs >/ki 

+rs < d^ijji, K-ijj^> /(/C2 - ki) if k2 = 3ki 

-ri < d^^i,K-^2 > /2/ci 



< d^ipi, K-h^> < dy'il^i, K-h4> T4 J 
h + h + ^<^^^i'-^^ > 

< dyipi, K ■ {hi - hi) > < dy'4)i, K ■ hi > 



ki 



ki - k2 



(44) 



(45) 



(46) 



Q3(o) 
Qi(o) 



< dyipi, K ■ {h; - hi) > , Tg 



ki 



ZK2 



<dJl,K-{hl-h2)> VI 



ko 



2k2 



where Vj = J2sI dvhf\v), and 



(47) 
(48) 



{^2,^^{r)) = Pi{0)B\A\'+p2{0)B\B\' + 5k„2kM0)A'' + 5kM{0)A' (49) 
\P2{0)AB' + {Qi{0) + P^{0))A*\B\' + Q^{0){A*)'B* + Qi{0)A*\A\''^ 



with coefficients 



Pi(0) = 



< dy'ii)2, K • (/i3 - ^5) > , < K-hi> Ta ~ 

H ; hi — < Oyijj2, n ■ ill > 



ki 



ko 



k2 - ki 



k2 + ki 



< 9^V'2, K- {h2- he) > +— < dy'il)2, k • V'2 > 



(50) 
(51) 
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-i < dvip2, K- ipi > /ki 



if ko 



2ki 



g(0) 



(52) 



-i < dyip2, K ■ >/ki + T4 < dyil)2, K- ipi > /2ki if k2 = 3ki. 
In (|4l|) and (^91) we abbreviate the notation of Section || letting rj{0) denote d^.r^O), -Pj(O) 
denote d^^P{0) and so forth. For an instabihty with complex eigenvalues and a reflection 
symmetry, the wavenumbers satisfy k2 = —ki, and the terms in r2(0) and J5i(0) are 
omitted. In addition, the reflection symmetry (p6D implies various identities: ri(0) = ^2(0), 
r2(0) = Pi(0), hi{v) = h2{-v), h-s^v) = hs^-v)*, hi{v) = ha^-v), and h5{v) = h^i-v), and 
these relate the cubic terms in ( ^T]) and P9| ) as shown. 

The coefficients hi{v) follow from the second-order terms in (pG]). On the left-hand side 
of fEBD, we have 



dS_ 
'dt 



dH : dH . 



dH ^ ^ dH ^ ^ 
-A,A + -A25 + cc 



+ C?(3), 



(53) 



with the partial derivatives evaluated using (^0|). On the right-hand side of (^6|), the leading 
terms are found in CH and 



^^^K ■ d^{i){ - ipi) + ^i^K ■ 9^(?/'2 - ^2) 
ki k2 

[ ki k2 \ ki k2 



(54) 



1 ) ^i{ki+k2)x 



V h k2 



i(fei —k2)x 



+ cc} + 0{3). 



The second-order solution of 



determines the coefficients in 



hiv) 



in ■ dy^ipl - ipi) 



h2{v) 



IK ■ (9„(^2 - ^2) 



271 fci 272 /C2 

iRk.^kAK + A2)/3(^^) h^{v) =iR2k,i2X,)h{v) 
iRk^+kA^i + A2)/5(^^) h(i{v) = iR2k2{'^X2)h{v) 



(55) 
(56) 
(57) 



where Ri{w) denotes the resolvent operator (|l5|) , and 
K ■ dyipl K ■ d^il)2 



hv) 



k2 



+ is{0) ipi6k2,2ki 



hiv) 



ko 



(58) 



14 



h{v) = — h I g(0) 42.2fci /el^^) = — ; ■ (59) 

These expressions are valid for each of the three instabihties we consider, with one exception. 
In the case of complex eigenvalues and reflection symmetry, when ki + k2 = 0, the result for 
/15(f) in ( pTl) is replaced by 

h{v) = (fci + k2 = 0) (60) 



in the notation of Section p.2.2 . 



Note that ( |55D implies Fi = r2 = in general, and (0) forces Fs = for this specific 
type of instability. Following (I), we can re-express (F^, hj), j = 3, 4, 5, 6 in (|55| ) - (0) more 
conveniently as 

where 

h = k2- ki, Z3 = k = '2ki, = Zi; (62) 

k2 - ki 

k = k2 + ki, Z5 = ^2^2 + l^i^i ■ ^ 2k2, zq = Z2. (63) 

K2 + «;i 

2.4 Previous results for the single mode instabilities 

For A = in (|30|) - ([3T|), we recover the single mode instability i? = p(|i?p)i? = [A2 + 
P2(0)|-Bp + ■ ■ \B previously considered JI], ||; the asymptotic form of the cubic coefficient is 

P2(0) = \\c{k2,Z2)-l2d{k2,Z2)^0{:il)\ (64) 

with the nonsingular functions c and d defined by 



c(.2,.2) - -iXT^E^Hl - .(^)^) Im (£ ..-^ 1 (65) 



-5]>)(i -/€«') imf r 



^fc2(^2)V" ' i-c^ (t; - 2:2) 
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where the primed species sum omits the electrons. |1|] For ni > 1 in pO| ) - (|3lD, setting 
-8 = determines a single mode instability at the longer wavelength: A = r{\A\'^)A = 
[Ai + ri(0)|74p + ■ ■ -jA. The asymptotic form of ri(0) follows from (|5^) - ( pBD with the 
replacements (72,^25-22) {li,ki,Zi). We also reduce to this previously studied case on 



setting either A = or i? = in equation (|38D ; these two limits correspond to single- 
mode instabilities in the form of travelling waves propagating in the positive and negative x 
directions, respectively. 

3 Singularities in the mode-mode couplings 

The central result of the single mode analysis proves that setting jSj = 5/2 in ([^) yields 
rescaled amplitude equations for a(r) and b{r) that are finite to all orders. When there 
are two unstable modes additional coupling terms between the modes are present that have 
not been previously considered. In particular in the expansions ( ^T]) and ( P9| ) we find four 



such couplings, (g(0), s(0),pi(0), r2(0)), for instabilities with k2 > ki > (cf. Sections |2.1| - 
T2A\) and five couplings, {pi{0),Qi{0), P2{0),Qs{0),Qi{0) +^3(0)), for an instability with 



^2 + ^1 = (cf. Section p.2.2|) . The singularities of such mode coupling terms will determine 
if the presence of a second unstable wave can alter the nonlinear scaling associated with 
single wave instabilities. 

This question needs to be formulated carefully to avoid "trivial" limits since there are 
now two distinct linear growth rates 71 and 72. If one growth rate vanishes while the sec- 
ond remains bounded above zero, then only the singularities associated with the resonant 
denominators of the first mode will emerge. This effectively recovers the singularity struc- 
ture of the one mode problem even though both mode amplitudes are non-zero. The more 
interesting limit arises when all resonant denominators come into play which requires both 
growth rates to vanish simultaneously. Thus we set 71 = 72 = 7 in the following discussion; 
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in practice this arrangement would be hard to reahze experimentally but could be achieved 
in a numerical simulation by simultaneously adjusting the parameters of the equilibrium and 
the length of the system. 

The origin of the singularities is the same as in the one mode problem: poles from resonant 
denominators can straddle the contour of integration and produce pinching singularities as 
7 — * 0+. Following the methodology of (I), the worst possible singularity of a given integral 
can be estimated by simply counting the total number of poles (including multiplicity). This 
gives an upper bound on the possible divergence of a given integral which must be checked 
by a detailed evaluation once the most divergent integrals have been identified. Ultimately 
we are most interested in possible singularities that are stronger than those already identified 
in the single-mode subsystem ([6^). For example, third order terms with divergences weaker 
than the 7"^ singularity of (|64D are clearly sub-dominant and cannot alter the (3j = 5/2 
scaling forced by the single mode singularities. 

In Sections P]T| - 0|, we assume the single mode problems exhibit the 7""^ cubic singularity 



in (H), i.e. 

^K(^)(l-K(^)')r^«(t;,)^0, J = 1,2 (67) 

s 

for modes ki and ^2, respectively. This is simply the condition that in (^) c{kj, Vj + iO) 7^ 0, 
j = 1,2. Special limits, such as infinitely massive ions, for which the single mode system 
is less singular are discussed in Section Tables 1 and 2 provide a summary of the 

asymptotic behavior of the second and third order coupling coefficients. 

3.1 The universal couplings pi{0) and r2(0) 

The coefficients r2(0) and pi{0) are present at third order for each instability, and their 
asymptotic singularities are fundamental. We discuss the three types of instabilities sepa- 
rately. 
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3.1.1 -Fo(f) 7^ Fo(— f), complex eigenvalues 



For instabilities without reflection symmetry (Section |2.1|) , we have ^2 + > for two 
positive and unequal wavenumbers. We first identify integrals in r2(0) and pi{0) with poles 
above and below the contour; terms without this feature are manifestly free of pinching 
singularities and will be finite as 7 ^ 0^. In addition, singularities weaker than 7"^ are 
sub- dominant. 

From (|62D - (|63D, lm{zj) > for j = 1, . . . , 6 so ipi{v), ip2{v) have poles only in the upper- 
half plane as do the conjugated adjoints ipj^v)*; thus < dvi'i, ■ ip2 > and < dv%l)2, k ■ ipi > 
are nonsingular. Similarly, /5(f) and h^{v) have poles only in the upper-half plane, thus 
the integral in (^) for Fs is nonsingular; in addition, IS^uX^^) = ^k2+kA^!3) is of order unity 
as 7 — since there are no modes at wavenumber k2 + ki by assumption. Hence F5 and 
< dv'ipi, /is > are both nonsingular. The integrals < d^ifji, k ■ > and < d^%lj2, n ■ ipt > are 
nonsingular unless Vi = f 2 in which case 7"^ is their worst possible divergence, so the terms 
involving F5 are sub- dominant. 

Similarly in 



< d^tpi, K- h2>-- 



27^2 



< dv^pi, ■ dyip2 > - < dyljji, ■ dvll)2 > 



(68) 



the second integral is nonsingular and the first integral, < d^tpi, ■ dy'il)2 >, is nonsingular 
unless f 2 = f 1 in which case there is a 7"^ singularity. This conclusion applies equally to the 
corresponding terms in < dvip2i i^-hi >; thus these terms contribute at most a 7"*^ singularity 
to r2(0) and Pi{0), respectively. 

Discarding these terms, we must still consider the asymptotic behavior due to and F3: 



r2(0) 
Pi(0) 



< dvipi, K-hl> 

< dy^2, K- h3> 



k2 - ki 



+ 



k2 - ki 



< dy'llj2, K,-ll>l> 



+ ■ 



+ 



(69) 
(70) 
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where 



1//3 \ r dvEs 



^hi^s) J J -00 V - Z3 



(71) 



r* 2 

<d,^uK-h;> = < a.^i, ' 3 > -r* < d.^Pu " ' \ > (72) 
<d,i>2,f^-h> = < d,i>2, ' , > -r3 < d,i)2, " "^'\ > . (73) 

h{v - Z3) [v - Z3) 



For (p9| ) - (|70D, we discuss the instabihties with ^2 7^ 2A;i and A;2 = ^ki separately; the latter 
case is more interesting since there can be new singularities when the modes have the same 
phase velocity, e.g. in a beam-plasma instability with a sufficiently cold beam. |T7| 



When k2 7^ 2ki, the third term in (fflf) is absent and the second term is manifestly 
nonsingular; from (|6^), 

^3 = ; + « 1 r ' (74) 



k2- ki \k2-k1, 

so the ffist term in (ffl]) has a pinching singularity only if f 2 = Vi and this possibility yields 
a divergence of 7"^. Also, A;3(2;3) = Kk,^_i^^{z^) is of order unity as 7 — > since there are 
no roots with Im(z) > for wavenumbers other than ^2 and ki\ hence the Fa terms in ( |69D 
- (^) are sub-dominant. The second term in ([73|) exhibits only the singularities of F3 and 
may be dropped, while the second term in ([7^ ) has an additional factor with a pinching 
singularity similar to F3, except that the roles of zi and z^ are reversed, and gives an overall 
divergence of at most 7"^ (when V2 = Vi). The remaining terms in ([7^) - ( |75D are 

h{v-zl) k2h{v-zl) kih{v-zl) 

< 77 r > = < d^il)2, j-r-, ^ > - < d.,i)2, J— >; (76) 

h[v - z^) k2k{v-Z3) kih{v - Z3) 

in all cases there are five poles in the integrand and hence a maximum possible 7"^ divergence. 
A closer examination shows that, if V2 = fi, this singularity is realized by the ffist term in 
([761) and both terms in (|75D; in any event we do not encounter a singularity in r2(0) or pi{0) 
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that dominates those found at third order in the single mode equations, i.e. a singularity 
stronger than 7"'^. 

When k2 = 2ki, Ai^iz^) = Ak^iz^) with 

zs = {2v2-v^)+^(^^y (77) 

Thus when V2 = vi we have Z3 = zi + i'-y/ki, and Ai^^zs) is 0(7) as 7 — 0"^; if V2 7^ vi, Ai^^z^) 
is still of order one. Thus with equal phase velocities, the singularity of T3 is increased to 
7~^ (including the presence of the third term in (|7l| ) which has a 7"^ divergence in s(0)); 



nevertheless the terms in (|69| ) - (|70D are still sub-dominant. The second term in ( [73D 
can be neglected for similar reasons. However, the second term in ([72|) now has an apparent 
divergence of 7"^ when V2 = Vi, such a singularity would not be absorbed by the scalings used 
to remove the 7"^ singularities characteristic of the single mode problem. We will evaluate 
this term more precisely shortly. The first terms in ([7^ ) - { \^^ ) contain extra contributions 
from the third term in (Wl\): 



<^.^P^,rT^^> = < d.^Pu ,f ; >-< d.^„ • > (78) 

h{v - z^) k2h{v - zl) kih{v - zl) 

h{v-zs) k2k{v - Z3) kil^^v-Z'i) 

h[v - 2:3) 



Our previous discussion of the first two terms in ([78| ) and (|79D is unchanged, and the new 
term in ([79| ) shows only the 7"^ singularity in s(0) (as determined from equation (^) below). 
This singularity is also present in the third term of (|78D, but now s(0)* multiplies an integral 
that has an apparent singularity of 7"^ giving a second term with overall 7"^ divergence. 

The foregoing discussion shows that there are new singularities in the coupling coefficients 
and that these singularities are most severe when ^2 = 2ki and vi = V2- In this resonant 
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case, the pi(0) singularity never exceeds ■y~'^, but we have identified two contributions to 
< dyipiyh^ > in the cubic coefficient r2(0) whose singularities appear to dominate the 7""^ 
divergence characteristic of the cubic terms of the single mode problem. We proceed to a 
detailed evaluation of these exceptional terms which shows that the 7"^ singularities are 
typically present. 

For the second term in ([7^) the calculation of the integrals using partial fraction expan- 
sions (cf. (I)) yields the following asymptotic forms when V2 = vi. 



1 



^'«:W(l + «:(^))r/g^(t;0 + O(7) 



MM) 



(80) 
(81) 



where the primed species summation excludes the electrons. In writing (|8l|) we have used 



the dispersion relation (piq) . Given the assumption (]67|) on the single mode problem, we 
expect Z]s''«*-'*'*(1 + i^^''^^)Vk'i i'^'i) 7^ to typically hold; hence the 7"^ singularity is generally 
present. Finally, the third term in (ffSD contains a 7"^ singularity from s(0)* (equation 
while a partial fraction expansion of the integral yields 



> 



^)E>ni + .^^)).i?(.0 + o(7) 



(82) 



{v - zi) r 

Thus the third term in (|78D also realizes a 7"^ singularity. These singularities require a shift 
in the scaling exponents that characterize the single mode instability; this point is discussed 
below in Section |[ 

3.1.2 -Fo(^) = real eigenvalues 



For reffection-symmetric instabilities with real eigenvalues (Section |2.2.1 ), we have Vi = 
V2 = and + ^1 > for two positive and unequal wavenumbers. The previous analysis 
is applicable here and we obtain the same conclusions with one modification. Since the 
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Table 1: Generic singularities of the couplings and scaling exponents 
Instability Resonance q s r\ r2 Pi pi {Pi, l32) 



Fo{v) ^ Foi-v) 


k2 


^ 2ki, V2 ^ vi 


- 






7" 


-4 


7" 


-1 


7" 


-1 


7" 


-4 


V 2 ' 2 


k2 > ki > 




V2 = f 1 








/ 


-4 


/ 


-4 


/ 


-4 


/ 


-4 


(5 5^ 
V2' 2/ 


Ai, A2 complex 
































k2 


= 2ki, V2 7^ Vi 


0{1) 




-2 


7' 


-4 


7' 


-2 


7' 


-2 


7' 


-4 


(5 5^ 
V 2 ' 2 ' 






V2 = Vi 




r 


-2 


7" 


-4 


7" 


-5 


7" 


-4 


7" 


-4 


(1,3) 




k2 


= 3ki, V2 ^ vi 


0{\) 


r 


-4 


7" 


-4 


7" 


-1 


7' 


-1 


7" 


-4 


f- -) 

V2' 2/ 






V2 = Vi 




r 


-4 


7" 


-4 


7" 


-4 


7" 


-4 


7" 


-4 


V2' 2/ 


Fo{v) = Foi-v) 




k2 ^ 2ki 








7" 


-4 


7" 


-4 


7" 


-4 


7" 


-4 


V2' 2^ 


k2> ki> 






























Ai, A2 real, 




k2 = 2ki 


0{\) 


r 


-2 


r 


-4 


r 


-5 


7" 


-4 


7" 


-4 


(1,3) 


multiplicity-two 


































k2 = Ski 


0(\) 


r 


-4 


7" 


-4 


7" 


-4 


7" 


-4 


7" 


-4 


f- -) 

V2' 2^ 



t = single mode 

condition vi = V2 is automatically satisfied, the "spatial resonance", k2 = 2ki, is sufficient 
to obtain the extra singularities noted above. These results are summarized in Table 1. 

3.1.3 Fo{v) = Fo{—v), complex eigenvalues 



For reflection-symmetric instabilities with complex eigenvalues (Section |2.2.2|) , we have Vi 



—V2 and k2 + ki = 0. These conditions rule out the presence of k2 = 2ki or k2 = 3ki 
resonances, and imply that pi{0) = ^2(0), where 



Pi(0) 



< (9i,V'2, ■ (^3 - ^5) > _^ < dyilj2, K- hi> ^ < d^tlj2, n-^i> 



(83) 



ki k2 k2 - ki 
with /15 defined in (|60|). 

From Zi = —Z2 = Zq we now have Zi in the upper-half plane (fci > 0) and Z2 in the 

lower-half plane (^2 < 0) along with = —i'-f/ki. The eigenfunctions ipiiv) and ipi{v)* 
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have poles in the upper-half plane that approach fi as 7 — > 0+, while the poles of 4'2{v) and 
il'2{'v)* are in the lower half-plane and approach f 2 as 7 — » 0+. The integrals < d^ipi, ■ V'2 > 
and < dvilj2, k ■ ipi > are still nonsingular although they now involve poles above and below 
the contour. The relation vi = —V2 ensures that no pinching singularity develops. Similarly, 
/5(f) and h^{v) have poles on each side of the contour but no pinch can develop in the 
integrals < d^ipi, k ■ h^> and < dyil)2, k ■ h^>. In 

< d^^Ji, K-h2>= [< ^^V"!, ■ dvip2 > - < <9^^i, 1^'^ ■ d^'ip2 >] , (84) 

the first integral < dyipi, ■ dytp2 > only has poles in the upper-half plane while the second 
integral has poles above and below the contour but no pinch; thus this term may also be 
dropped as well as the corresponding term < dyilJ2, n ■ hi > in pi{0). 

Discarding these manifestly nonsingular terms, we must still reconsider the remaining 



terms in (|69D - (^). The prior analysis given in Section 3.1.1 requires relatively minor 
modifications to allow for the relation vi = —V2 and the shifted location of the poles Z2 
and Z3. In (|7T]), A/3 (2:3) = A_2fci(-23) is of order one as 7 — *• 0"*" and the integrals are free 
of pinching singularities since Z3 is pure imaginary and can never converge to either zi or 
Z2. Thus the second terms in (^) - (|70D are nonsingular as are the second terms in (^) - 
([73|). The final integrals in the first terms of (|72D - ([73l), written out in (^) - (|76D , are also 
manifestly free of pinching singularities. In summary, the coupling coefficients Pi(0), r2(0) 
for instabilities with complex coefficients and reflection-symmetry are identical, and exhibit 
singularities due to the explicit 7"^ factor in < d^ipi, k ■ h2 > only. 

3.2 The couplings Qi(0), P2(0), QsiO) and Qi(0) + P^iO) 



For reflection-symmetric instabilities with complex eigenvalues (Section |2. 2. 21) , there are four 



additional 0(2) symmetric couplings at third order. Although these terms can be removed 
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Table 2: Generic singularities of the couplings and scaling exponents 
Instability Qi P2 Q3 (Qi + A) A Pi pi {f3u (^2] 



Fo{v) = Fo{-v) 0{1) 7-' 

k2 + ki = 
Ai complex, 
multiplicity-two 



7" 



^4 7"'" T"'^ 



'5 5^ 
^2' 2' 



t = single mode 



by a coordinate transformation to obtain the normal form in (pSf), it is important to consider 
their asymptotic behavior. 

The singularities for the couplings in this instablity are summarized in Table 2. ^2(0) 
is nonsingular and the singularities of Qi{0) and Q^iO) are sub-dominant. The strongest 
singularity occurs in Qi{0) + PsiO) due the integral 



< dvipi, K ■ hi >-- 



< dytpi, ■ dyipl > — < dv^i, ■ dvipi > 



(85) 



whose first term < dyipi, k"^ ■ dyipi > has a 7 singularity giving an overall singularity of 7 ^ 
for gi(0) + P3(0). 

3.3 The spatial resonances: q{0) and s(0) 

When k2 = 2ki and k2 = 3ki there are additional couplings at second and third order, 
respectively. We first consider g(0) and s(0) for the k2 = 2ki resonance, noting that the 
two integrals, < d^ipi, n ■ ip2 > and < dvip2, «^ ■ V^i >, have poles at zi and Z2 and are free 
of pinching singularities in all cases. Thus g(0) is nonsingular as is the second term in s(0). 
The first term in s(0) was evaluated in (I), 



2i 



KM) 



„ J-00 V — Zi 



16) 
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where the primed species sum omits the electrons; this determines the singularity of s(0), 



"^0) = ^2 



Tiki 



(87) 



For the ^2 = 3fci resonance g(0) is again readily seen to be nonsingular and we omit the 



details. For s(0) in (^J) the integral < d^'ipi, n ■ h\> contains a contribution, 



V — Z^ 



if V2 7^ Vi. When V2 = fi there are additional singularities, but none stronger than 7 ^. 
These conclusions are summarized in Table 1. 

3.4 Special limits: coupling singularities with fixed ions 



In the various exphcit asymptotic formulas, such as (0), (F0|) - (|82D , and (|87|) , the leading 



term vanishes if the ion masses are treated as infinite since n^^^ ^ in this limit. In (I), this 
suppression of the leading singularity was shown to be a general feature of the integrals that 
appearing in the amplitude equations. For the single mode instability, the cubic coefficient 
Pi(0) has a 7"^ singularity when the ions are fixed, and the modified single-mode scaling 
A(t) = 7^a(7t) suffices to render the amplitude expansion finite. ^ 

It is straightforward to adapt the results of the previous sections to the case of infinitely 
massive ions: with only a few exceptions among terms that are already sub-dominant the 
generic divergence is reduced by one factor of 7"^. For our purposes it is suffices to summarize 
in Table 3 the resulting changes to Table 1 when the electrons are the only mobile species. 

4 Nonlinear scalings 

In this section we make use of the leading order behavior of the coupling coefficients identified 
in the preceding section to determine the scaling of the saturation amplitude of the two 
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Table 3: Summary of the singularities of the couplings for fixed ions 



Instability Resonance q s rj r2 pi pi {Pi, P2) 



Fo{v) ^ Foi-v) 


k2 


^ 2ki, V2 ^ vi 


- 






7" 


-3 


7" 


-1 


7" 


-1 


7' 


-3 


(2,2) 


k2 > ki > 




f 2 = Vi 








/ 


-3 


/ 


-3 


/ 


-3 


/ 


-3 


(2,2) 


Ai, A2 complex 
































/C2 


= 2ki, V2 7^ Vi 


0{l) 




-1 




-3 


7~ 


-1 




^1 


7' 


-3 


(2,2) 






V2 = Vi 


0{l) 


r 


-1 


r 


-3 


7" 


-4 


7" 


-3 


7" 


-3 


(2,1) 




k2 


= 3ki, V2 ^ vi 


0{l) 


r 


-3 


r 


-3 


7" 


-1 


7" 


'1 


7' 


-3 


(2,2) 






V2 = Vi 


0{l) 


r 


-3 


7" 


-3 


7" 


-3 


7" 




7' 


-3 


(2,2) 


Foiv) = Foi-v) 




k2 ^ 2ki 








7" 


-3 


7" 


-3 


7" 


-3 


7" 


-3 


(2,2) 


k2> ki> 






























Ai, A2 real, 




k2 = 2ki 


0{l) 


r 


-1 


r 


-3 


r 


-4 


7" 


-3 


r 


-3 


(2,1) 


multiplicity-two 


































^2 = 3fci 


0{l) 


r 


-3 


7" 


-3 


7" 


-3 


7" 


-3 


r 


-3 


(2,2) 



t = single mode 



competing modes with the growth rate 7. We first consider this question for the coefficient 
singularities listed in Table 1, i.e. for instabilities with k2 > ki > 0. Following (I), we 
introduce scaled amplitudes A(t) = 7^^0(7^) exp(— ■j^^i(t)) and B(t) = ■y^^a{'yt) exp{—i92(t)) 
for 7 > and rewrite (p^ - (^) using the expansions in (^) and (^Of ): 

^ = a + 'j^'-^Re[s{0)e'^'''-''^]abSk„2k,+7"''-'Mrmy (89) 
ar 

+7'^^-^ Re[r2(0)]a + Re[s(0) e'^^''-''^]aH6k,,3k^ + ■■■ 

^ = 6 + 7'^^-^^-'Re[g(0)e-(2^-^^)]a2 4„2.,+7'^^-'Re[pi(0)]a2 6 (90) 
ar 

+^3/3i-A-i j^e[g(o) e-'^^''-'^'>]a^ 6k,,3k, + 7'^'"' Re[p2(0)] + ■ ■ ■ 
^ = c.i-7^^Im[s(0)e^(2^-^^)]64..2..-7'^^Im[ri(0)]a2 (91) 
-72/32 Im[r2(0)] 6=^ - 7^2+^^ lm[s(0) e'^^''~'^^] ab6k,,3k, + ■■■ 
26 



^ = ^2 - l''^-'^ lm[g(0) e-^(^^-^^)]^ - 7'"^ Mp^V (92) 

-7'^^^ Im[p2(0)] - r'^^-^' lm[g(0) e-(3e-^^)]^ 4„3fc, + ■ ■ ■ • 

If possible, the choice of /3i and /32 should be made so that each term is finite as 7 ^ 0+ and 
there is a formal balance between linear and nonlinear terms in (89) - (9C). 



In Table 1, we focus initially on the instabilities with 7^ Ik^. The 7"^ singularities of 
the single mode coefficients ri(0) and ^2(0) in (|89D - ( pOD require (5\ > 5/2 and P2 > 5/2. 
These exponents are large enough to ensure a finite limit for each term (^) - (P^; in fact 
the phase equations reduce to 9j = ojj + 0(7), j = 1,2. In (^) - (^0]), the minimal choice 
13 j = 5/2 suffices to control the singularities in the mode couplings r2(0) and pi(0) as well 
and the amplitude equations reduce to 

^ = a + Re[c(A;i, vi + tO)]a' + 7^ Re[r2(0)]a + 7" Re[s(0) e'^^''''''>]a^ b 6k„3k, + ■ ■ (93) 
ar 

^ = 6 + 7^ Re[pi(0)]a2 b + Re[c(fc2, t;2 + ^0)] b^ + 7^ Re[g(0) e-'^^''-'^^]a' 6k,,3k^ + ■ ■ (94) 



where c{kj,Zj) is defined in ([65|). In these variables the linear growth rates are unity and 
the single mode terms, and b^, are of 0{1) in 7. If V2 = vi, then the coupling coefficients 
7^^2(0) and 7^1(0) S'i's also of 0{1), but otherwise these terms are sub-dominant. When it 
is present, the resonant term, 7^s(0) e*^^^^"^^)^ jg formally of 0{1); however the phases in the 
exponential evolve on the fast time scale t = r/'y and therefore the exponential oscillates 
very rapidly as 7 — O"*" unless the linear frequencies are resonant: 3uJi = 002- Such a rapid 
oscillation allows the term to be time-averaged and neglected, but when "iuoi = uj2 the phase 
is stationary and this 0{1) term cannot be discarded by time- averaging. Note that, given the 
spatial resonance ^2 = 3^1, the frequency resonance is equivalent to assuming that the linear 
phase velocities are equal V2 = fi- The second resonant term, 7^g(0) e"**-'^^^"^^'', is formally 
of 0(7^) but nevertheless has a qualitatively important effect on the dynamics near 6 = as 
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discussed below. Thus the asymptotic equations (pSD - (p^ ) are correct when V2 = vi, but if 
f 2 7^ f 1 then their form simphfies to 

^ = a + Re[c(A;i,t;i + iO)]a^ H (95) 

ar 

— = 6 + Re[c(fc2,t;2 + ^0)]63 + yRe[g(0)e-^(3^^^-^2)]a3 4^,3,^ + ...; (96) 
ar 

this system describes the evolution of the waves as if they were non-interacting, indeed 
their mutual dynamics is essentially the "superposition" of the two single mode amplitude 
equations save for the very weak resonant term in (|96|) . For k2 7^ 2A;i, our principal conclusion 
is that no change from the scaling exponents predicted by the single mode singularities is 
indicated. 

The resonance ^2 = 2A;i is special: the s(0) resonant term now occurs at second order with 
a singularity and, if V2 = Vi, the singularity in r2(0) jumps to 7"^. In addition, t>2 = Vi 
implies a frequency resonance 2uji = 002 that suppresses the oscillation of the quadratic term. 
If we keep /?2 = 5/2 in (p9D, then this quadratic term retains a 7^^/^ singularity that cannot 
be erased by time averaging. Similarly the cubic term 7^ Re[r2(0)]a 6^ would retain a 7"^ 
singularity. In this situation, the scaling exponent for the short wavelength mode must be 
increased to /52 = 3 to obtain a sensible asymptotic limit. Now the amplitude equations 
reduce to 

^ = a + 72Re[s(0)e^(2ei-^2)]a6 + Re[c(fci,^;i + zO)]a3 + 7^Re[r2(0)]afe' + --- (97) 
ar 

^ = 6 + 7Re[g(0)e-'(2'i^'^)]a2 + 7^Re[j9i(0)]a26 + ---; (98) 
ar 

in these variables the growth rates are again unity and all terms shown are 0{1) in 7 
except for 7 Re[g(0) exp — i(2^^i — ^2)]^^^ which is 0(7) but is nevertheless important for the 
asymptotic dynamics as explained below. When ^2 = 2A;i but V2 7^ fi, the quadratic term in 



(p7|) can be removed by time- averaging and the singularity of r2(0) and pi(0) drops to 7 ^; 
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in this case the single mode scahngs suffice to control the singularities in ( ]89| ) - ( ]90| ) and the 
amplitude equations simplify to 

^ = a + Re[c{ki,vi + iO)]a^ + ■■■ (99) 
dr 

^ = b + ^Re[q{0)e-'^^''~''^]a^ + Re[c{k2,V2 + tO)]b^ + ■■■. (100) 
dr 

In the resonant regime, k2 = 2ki and V2 = Vi, the long wavelength mode saturates at 
a 0(7^/^) amplitude while the short wavelength mode saturates at a much smaller 0(7^) 
amplitude, i.e. the presence of the long mode results in a dramatic suppression of the 
amplitude of the short mode. Because of the (32 = 3 scaling required for the short wave 
mode the equation for 6(r) is dominated by the interaction with the long wave mode a(r). 
In particular in the rescaled equation (|98|) the single mode terms appear at higher order 
and hence are omitted. On the other hand we have retained the 0(7) term a^. Although 
formally small this term exerts a profound effect on the resulting dynamics. This is because 
it destroys the invariance of the subspace 6 = 0. This term therefore plays the role of a 
"noise" term that continuously forces the short mode. This observation holds equally for the 
resonant term in (0). 

The singularities summarized in Table 2 refer to an instability of a reflection-symmetric 
equilibrium involving complex eigenvalues. All the couplings in Table 2 are third order 
terms and none of the divergences are stronger than the 7"^ singularities of the single mode 
couplings. Thus the single mode exponents f3j = 5/2 will control all singularities in the 
amplitude equations to third order, and the rescaled equations have the form 

^ = a + Re[c{kuVi + iO)]a^ + -f^Re[e'^''+'^\Qi{0) + P3{0))]ba^ + ■■■ (101) 
dr 

^ = b + Re[c{k2,V2 + tO)]b^ + 7^Re[e'^''^+'-\QiiO) + PsiO))]ab^ + ---. (102) 
dr 

The joint conditions k2 + ki = and V2 = — f 1 characteristic of this instability imply u)2 = 1^1; 
hence the exponentials in the mode coupling terms are rapidly oscillating and we expect these 
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terms to have a negligible effect on the time- averaged evolution. In effect, the coupled system 
( |101| ) - ( |102| ) reduces to ( ]95| ) - (0) and we recover the superposition picture of the two single 
mode amplitude equations but without the resonant a? term. 

If the ions are taken to be infinitely massive, then Table 3 replaces Table 2 and the 
scaling associated with the single mode instabilities becomes (3j = 2. The re-analysis of 



(p9|) - (|9(]|) leads similar conclusions although with quantitative differences. Only for the 
resonance k2 = 2ki and V2 = fi are the single mode scalings altered and again the short 
wavelength mode is suppressed; in this case we find P2 increases to P2 = 5/2. When V2 7^ vi 



we again find that the amplitude equations decouple as in ( p5D - (p6|) . 

5 Discussion 

In this paper we have investigated the interaction between two weakly unstable electrostatic 
waves in an unmagnetized plasma. Such growing modes are defined unambiguously by 
eigenfunctions of the linear theory; Landau-damped collective modes are properly thought 
of as part of the continuum and hence are not decaying "modes" in the present sense. 0] 
For the growing modes discussed here, the theory provides a self-consistent description that 
includes the nonlinear effects of wave-particle resonance. Such a formulation, while more 
complex, captures phenomena that are absent from conventional formulations involving wave 



interactions. 11 



The wave-wave interaction is described by coupled amplitude equations consistent with 
the assumed translation invariance of the system and any symmetry of the unstable velocity 
distribution function. The coupling coefficients exhibit new singularities in the weakly unsta- 
ble limit, but except in the case of the two-to-one spatial resonance do not alter the scaling 
for the saturation amplitude identified in single mode theory, at least through third order 
in the mode amplitudes. In the special but important case of two-to-one spatial resonance 
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the overlap of velocity resonances and spatial resonance modifies the scaling, resulting in a 
dramatic suppression of the short wave mode. 

In the absence of velocity overlap and the two-to-one spatial resonance the single mode 
scaling shows that the two waves evolve as a simple superposition of the individual instabil- 
ities. This picture is qualitatively consistent with the numerical results obtained by Demeio 



and Zweifel for beam-plasma instabilities with reflection symmetry. [|10[ In addition, it agrees 
more quantitatively with the analysis of Buchanan and Doming who constructed superpo- 
sitions of BGK modes as candidates for the asymptotic states produced in the numerical 



simulations. |]T2[ In particular, these authors found that a consistent construction to leading 
order in the amplitude of the individual BGK modes required unequal phase velocities; an 
assumption analogous to our V2 7^ f 1 condition. Our results suggest a precise connection 
between the initial value problem for the unstable waves and the superposed BGK states of 
Buchanan and Doming. 

The significance of resonance overlap, such as V2 = f 1, in the single particle phase space is 
well-established from studies of particle motion in fields produced by large amplitude waves. 
In these studies the appearance of chaotic particle trajectories is investigated, but the self- 



consistent modification of the wave evolution by the particles is routinely neglected. [|13 



14| , [iq] The situation we consider does not allow this simplifying approximation since the 



resonant particles drive the linear instability and also dominate the nonlinear evolution of 
the waves. It is striking to discover that resonance ovelap has a profound effect on the 
dynamics of the waves, in addition to its better known consequences for the associated 
particle dynamics. 

Our prediction that resonant interaction with a longer wavelength mode can modify 
the nonlinear scaling of a short wavelength mode may be amenable to experimental test. 
The single mode scaling for fixed ions was detected experimentally in measurements on an 
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electron beam injected through a travelhng wave tube. [0] The tube plays the role of the non- 
resonance electrons and supports a propagating wave which couples to the resonant particles 
in the beam. If the electron beam is sufficiently cold the unstable waves approximately 
satisfy the approximate dispersion relation uj{k) = Vpk with constant Vp and hence always 
have equal phase velocities.] 17] Under these conditions, one should measure the scaling of a 



single mode launched at k2 and then repeat the measurement when a second wave is launched 
simultaneously at ki = k2/2. 

The possibility remains that singularities in the higher order coupling coefficients could 
modify these conclusions and force new scalings for other spatial resonances as well. How- 
ever, in the study of singular amplitude equations in other problems it is commonly found 
that the dominant singularities appear in the low order nonlinear terms. 1^ In addition, 
the fact that the k2 = 3ki spatial resonance does not shift the scalings may signify that 
resonances other than k2 = 2ki will generally leave the single mode scalings undisturbed. 
These speculations can be tested by examining the singularities in the amplitude expansions 
to all orders. Such an investigation may be feasible using the techniques of (I). 
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